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\0 ■ Abstract 

Critical phenomena and universality behavior of ferromagnetic thin films described by a spin- 

' (— l '' 
Q \ 1 Blume-Capel Hamiltonian has been examined for various thickness values ranging from 3 to 

(D ■ 

^ ' 40 layers. Using effective field theory, we have found that crystal field interactions significantly 



affects the critical value of surface to bulk ratio of exchange interactions Re at which the critical 
"^ \ temperature becomes independent of film thickness L. Moreover, we have extracted the shift 

exponent A from computed data. Based on the results, we have shown that in the presence of surface 
— L ■ exchange enhancement, the system may exhibit a dimensional crossover. We have also found that 

Q . presence of crystal field interactions does not affect the value of A. Hence, a ferromagnetic spin-1/2 



thin film is in the same universality class with its spin-1 counterpart. 
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I. INTRODUCTION 



In recent years, influences of surface effects on the magnetic properties of finite systems 



QE 



sucli as ferromagnetic tliin films liave attracted a considerable amount of interest [1|, |2| . In 
addition to remarkable theoretical efforts, preparation of thin films by depositing a magnetic 
material on a non-magnetic substrate became experimentally accessible even in the mono- 
layer limit with the development of modern vacuum techniques including molecular beam 
epitaxy. Due to the presence of free surfaces, magnetic properties of thin films may differ 
from those of bulk materials. This difference mainly originates from a number of physical 
phenomena. Namely, the surface atoms have a lower symmetry in comparison with that of 
the inner atoms [3], and the exchange interactions between the surface atoms may be dif- 
ferent from those between the corresponding bulk counterparts. As a consequence of these 
facts, the surface may exhibit an ordered phase even if the bulk itself is disordered which 



has already been experimentally observed 



4- 



In this context, an extraordinary case is defined as the transition at which the surface 
becomes disordered at a particular temperature T^ which is larger than the bulk transition 
temperature TJ!. From the academic point of view, due to the fact that many thin films 







such as the Fe/Ag(100) system [7| exhibit a strong uniaxial anisotropy, phase transition 
characteristics of thin ferromagnetic films are often modeled by several extensions of an 
Ising type spin Hamiltonian [8]. It is theoretically predicted that there exists a critical 
value of surface to bulk ratio of exchange interactions Re above which the surface effects 
are dominant and the transition temperature of the entire film is determined by the surface 
magnetization whereas below Re, the transition characteristics of the film are governed by 
the bulk magnetization. The critical value Re itself is called as the special point, and the 
numerical value of this point has been examined within various theoretical techniques for 



spin-1/2 case 9l-|l6l|. Among these works, within the framework of effective field theory 



(EFT), Sarmento et al. 1^ clarified that a transverse field in the surface layer causes the 
critical value of the surface exchange enhancement Re to move to a higher value whereas 
the presence of a bulk transverse field causes Re to decrease to a lower value. 



The problem has also been handled for higher spins using a number of techniques 17H24 1 . 



For instance, the phase transition properties of a diluted Blume-Capel ferromagnetic film 
with S" = 1 in a transverse field have been investigated by EFT [l8|. Using the same method. 



phase diagrams, the layer longitudinal magnetizations and quadrupolar moments of a spin-1 
film have been examined as functions of the ratio of the surface exchange interactions to 
the bulk ones, transverse fields, and film thickness 19|, |22|. In addition, using extensive 



Monte Carlo (MC) simulations, the effect of surface exchange enhancement on ultrathin 



spin-1 films has been studied by Tucker 



20|, and it was concluded that the R^. value is spin 



dependent. However, in a recent work 23[, using MC simulations, the influence of crystal- 
field interaction (or single ion anisotropy) on the critical behavior of a magnetic spin-1 
film has been studied, and it has been argued that Re is independent of the crystal-field 

n 

interaction. It is clear that this latter result apparently conflicts with [20], and according 
to us it deserves particular attention. Moreover, we see from above discussions that the 
theoretical studies based on EFT are completely focused on the effect of surface and bulk 
transverse fields in the absence of crystal fields. 

On the other hand, theoretical and experimental investigations are also focused on the 
finite size shift of the critical temperature of the film as a function of its thickness which is 
characterized by a shift exponent A. A number of experimental studies have been devoted 



to determine the value of A for various thin fi 
shift exponent extends from ~ 1 to 3.15 |25 



m samples, and it has been concluded that the 
428l|. Since the exponent A is directly related 



on the bulk correlation length exponent as v^ = 1/A [29\, within the accuracy of Ising-type 
models, it can be mentioned that a sample of thin film for which the exponent A is close to 
unity exhibits a two dimensional character whereas as the value of the exponent becomes 
larger than unity then the system shows a three dimensional character. Theoretically, the 
exponent A has been extracted for some certain models with a wide variety of techniques. For 



instance, using the high temperature series expansion 



that the estimated value of A for ferromagnetic Ising 30 



: TSE) method, it has been shown 



3l| and Heisenberg 32] thin films 



severely depends on whether a periodic or free boundary condition was considered in the 



surface. This result has also been verified within the resolution of MC simulations [331-136]. 
Moreover, the value of the extracted exponent is also very sensitive to the lattice geometry 
37]. 

Under certain circumstances, universal behavior of a thin film system may experience a 
dimensional crossover. Such a phenomenon has been experimentally observed as the film 
thickness is varied in ultrathin Ni(lll) films on W(llO) 38|], and epitaxial thin films of Co, 
Ni, and their alloys grown on Cu(lOO) and Cu(lll) 39|]. Previous MC simulations lOj also 



predict that the exponent A may vary continuously with surface exchange Js in the range 
R < Re which also indicates the occurrence of a dimensional crossover between the surface 
value and bulk value. 

As seen in the above discussions, there are some important questions which remain to 
be answered regarding the phase transition characteristics of thin ferromagnetic films. For 
example, the controversy on the crystal field dependence of special point Re should be 
resolved. In addition, it is not clear whether the presence of a crystal field interaction in 
the model affects the universality behavior of the system or not. In order to clarify these 
issues, we intent to study the spin-1 ferromagnetic thin film model in the presence of crystal 
field interactions within the EFT framework |40]. Despite its mathematical simplicity, EFT 
systematically includes the single-site correlations in the calculations, hence the obtained 
results are expected to be more accurate than those obtained by conventional mean field 
theory (MET). 

The organization of the paper is as follows: In Sec. [TTlwe briefly present the formulations. 
The results and discussions are presented in Sec. IIIIl and finally Sec. HV] contains our 
conclusions. 

II. FORMULATION 

We consider a ferromagnetic thin film with thickness L described by conventional Blume- 



Capel Hamiltonian 41 1 



n = -J2J^AS,-DJ2is,)'', (1) 

where Jij = Jg if the lattice sites i and j belong to one of the two surfaces of the film, 
otherwise we have Jij = Jb where Jg and Jb denote the ferromagnetic surface and bulk 
exchange interactions, respectively. The first term in Eq. ([1]) is a summation over the 
nearest-neighbor spins with Si = ±1,0 and the second term represents the crystal field 
energy on the lattice. 

The magnetizations and quadrupolar moments (i.e. rrii = {Si) and qi = ((Sj)^),? = 
1, ...,L) perpendicular to the surface of the film corresponding to L parallel distinct layers 
can be obtained by conventional EFT formulation based on differential operator technique 



and decoupling approximation (DA) J42| 



mi 
m„ 



= [1 + miBi + qi{Ai - \)X\\ + m^B^ + ^2(^2 - l)]Fi(x)U=o, 
= [1 + mp52 + gp(A2 - 1)]"[1 + mp_i52 + q^-x{A^ - 1)], 
+ [1 + mp+iB2 + qp+i{A2 - l)]Fi{x)\^=o, 
rriL = [l + mLBi + qL{Ai-l)Y[l + mL^iB2 + qL^i{A2-l)]Fi{x)\^=o, 
qi = [l + miBi + qi{Ai-l)Y[l+m2B2 + q2{A2-l)]F2{x)\^=o, 
qp = [1 + mpB2 + qp{A2 - 1)]"[1 + mp-i^a + 9^-1(^2 - 1)], 

+ [1 + mp+iB2 + qp+i{A2 - l)]F2(x)|a;=o, 
qL = [I + ttilBi + qL{Ai ~ 1)Y[1 + mL_iB2 + qL-iiA2 - l)]F2{x)\,=o, 

(2) 

where 2<p<L — l,zis the coordination number of the lattice, and the coefficients Ai and 
Bi are defined as Ai = cosh(JsV), A2 = cosh(JftV), Bi = sinh(JsV) and B2 = sinh(JftV). In 
the present work, we will focus on the ferromagnetic films in a simple cubic lattice structure, 
i.e. z = 6. The functions Fi{x) and F2{x) in Eq. ([2]) are then given by 

2 sinh(/3x) 



Fi{x) 
F2{x) 



2cosh{f3x)+exp{-f3D)' 
2 cosh(/3x) 



2cosh(/3x)+exp(-/3L))^ 
where f3 is the inverse of the reduced temperature. 



With the help of the Binomial expansion, Eq. ([2]) can be written as follows: 

z i 1 k 

i=0 j=0 fc=0 1=0 
z i 1 k 1 n 

i=0 j=0 fc=0 «=0 n=0 t=0 

i=0 j=0 fc=0 «=0 
z « 1 fc 

i=0 j=0 k=0 1=0 
z i 1 k 1 n 

i=0 j=0 fe=0 i=0 n=0 i=0 
i=0 jr=0 fe=0 «=0 



with the coefficients 



(3) 



V / V / "^"^ ^"° y 3; y y y 

x(-l)^+^+"-^-'-*-^e,(0,0,x,j + / + t), (4) 

where 

e„(A;,m,/,n) = AtA^5lfi2"F«(a^)U=o, « = 1,2. (5) 

Consequently, applying the Binomial expansion in Eq. ([5]) yields 

r=0 s=0 t=0 i>=0 y '^ y V '^ / \ ^ J V ^ / 

X exp[(2r + 2s - A; - /) J^V] exp[(2t + 2v - m - n)JbV]Fa{x)\^=o. (6) 

Using Eq. (Q, the coefficients in Eq. (^ can be numerically evaluated with the help of 
the relation exp(aV)/(x) = f{x + a) for an arbitrary a. Once the coefficients in Eq. ([6]) 

6 



are numerically evaluated, we obtain a system of coupled non-linear equations from Eq. ([3]) 
which contains 2L unknowns. The longitudinal magnetization rrii, as well as the quadrupolar 
moment Qi of each layer can be obtained from numerical solution of Eq. (^. Then, the total 
magnetization and quadrupolar moment of the entire system can be defined as 

1 ^ 1 ' 

m = -^mi, q=-^qi. (7) 

Since, the magnetization of the entire system is close to zero in the vicinity of the second 
order phase transition, the transition temperature can be obtained by linearizing Eq. ([3]), 
i.e. 

2 1 Z 

i=l fc=0 i=0 

z 1 z 1 1 



nip 



Y, Yl ^2'^ (^' 0' 1' 1' ^' 0)"^P-i«+i + J2J2Y1 ^2" ih 1, k, 0, n, 0)mpq'^_,q;-\;_,, 

i=0 n=0 i=l fc=0 ra=0 

z 1 

+ 5Z 5^ ^^'H^, 0, k, 0, 1, l)mp+igj_ig; 



i=0 A:=0 
z 1 2 

j=l fc=0 i=0 

2 1 



^1 = EE^i'^(^'0'^' 0)^1^2, 

i=0 fc=0 

2 1 1 

^p = EEE^2'^(^'0'^'0'^'0)^J-i^x+i 

j=0 fc=0 n=0 
z 1 

i=0 fe=0 

Critical temperature as a function of the system parameters can be determined from 
det(yl) = where A is the coefficients matrix of the set of 2L linear equations in Eq. (|8]). 
We note that the determination of the transition temperature should be treated carefully 
since as it was previously stated in Ref. IJ], from the many formal solutions of det(y4) = 0, 
we have to choose the one corresponding to the highest possible transition temperature. 

According to the finite-size scaling theory |29| , the deviation of the thickness dependent 
critical temperature Tc{L) of a thin ferromagnetic film from the bulk critical temperature 
Tc(oo) can be measured in terms of a scaling relation 

e = l- T,(L)/T,(oo) oc L-\ (9) 



for sufficiently thicker ffims where A is called the shift exponent which is related to the 
correlation length exponent of the bulk system as A = 1/vb- The exponent A can be extracted 
from numerical data by plotting e versus L curves for sufficiently thick films in a log — log 
scale then fitting the resultant curve using the standard linear regression method. 

III. RESULTS AND DISCUSSION 



In this section, we will discuss the effect of the presence of a surface, as well as the single- 
ion anisotropy on the universality features of layering transition characteristics of the system. 
At this point, we note that the value of the bulk exchange interaction J^ is fixed to unity, 
and we also use the normalized surface to bulk ratio of exchange interactions R = Js/Jb, 
as well as the reduced single-ion anisotropy energy A = D/Jf, and reduced temperature 
ksT / Jb throughout the calculations. For simplicity, the exchange couplings are restricted 
to the ferromagnetic case. 
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FIG. 1: (Color online) (a) Phase diagrams of Blume-Capel thin film in a {ksTc/ Jb — R) plane 
corresponding to the highly anisotropic limit A ^^ oo for various film thickness L. (b) Variation of 
the shift exponent A with surface to bulk ratio of exchange couplings R corresponding to A ^^ oo 
limit. 



In order to provide a testing ground for our calculations, we have primarily studied the 
phase diagrams in a (kBTc/Jb — R) plane corresponding to a highly anisotropic case A — )• oo 
with various film thickness L in Fig. [U^a). It is a well known fact that the behavior of 



the system reduces to that of a spin-1/2 thin film in this hmit. Namely, the special value 
of the surface to bulk ratio of exchange interactions at which the critical temperature is 
independent of thickness L is found as Re = 1.3068, and the corresponding bulk transition 
temperature is kfiTr/Jb = 5.0732. These results are identical to those obtained in previous 



works using EFT [ll|, llJ]. It is well established that for R < R^, the bulk dominates against 
the surface whereas for R > Re, the surface effects become prominent. In order to examine 
the universality behavior in this limiting case (i.e. A — )■ oo), we have extracted the shift 
exponent A by computing the transition temperature ksTc/Jb as a function of thickness L 
according to Eq. (jH]). In order to precisely cover the critical reg ion, the obtained data have 
been fitted for those providing the condition 0.01 < £ < 0.1 [33| which requires to consider 
the transition temperatures of the films for L > 10 in fitting procedure. The results are 
depicted in Fig. [^b). For R = 0.5, at which the bulk highly dominates against the surface, 
the estimated value of the shift exponent is A = 2.01 indicating the fact that the system 
exhibits likely a three dimensional character. On the other hand, in the vicinity of the special 
point R = 1.3 ~ Re, the surface effects become dominant, and we obtain A = 1.08 which 
reveals a two dimensional character even for thicker films. We can see from this result that 
as R increases continuously from zero to Re, the system may exhibit a dimensional crossover 
due to the presence of surface exchange enhancement. In addition, the value A = 1.87 
corresponding to i? = 1 can be compared with the results A = 1.56 and A = 2.002 obtained 
by MC simulations [33[ and MET |36|, respectively. 

In Ref. [20|, surface exchange enhancement effects on ultrathin spin-1 Ising films have 
been examined, and based on the results obtained within MC simulations, it has been 
concluded that Re is a spin dependent parameter which has a value lower than that obtained 



for a spin-1/2 system 13|. On the contrary, Ref. 23|] suggests that Re was independent 
of crystal field interactions which gives rise to a controversial situation. To the best of our 
knowledge, this case has not yet been examined in the literature. Therefore, in Fig. [2l we 
depict the phase diagrams in a (/cbTc/ Jf, — R) plane with various thickness L corresponding 
to some selected values of crystal field interactions A. It is clear from Fig. |2] that Re has 
its minimum value for a sufficiently negative crystal field value. For A = 0.0, we have 
Re = 1.2932 which can be compared with Re = 1.45 of MC simulations [20] whereas MET 

predicts a spin independent value. 

Dependencies of the special point Re and corresponding bulk transition temperature 
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FIG. 2: (Color online) Phase diagrams of Blume-Capel thin film in a {ksTc/ Jh — R) plane for 
various film thickness L corresponding to (a) A = —2.0, (b) A = —1.5, (c) A = —1.0, and (d) 
A = 0.0 

ksTc/Jb on crystal field interaction A have been clarified and presented in Fig. [3l As seen 
from Figs. EJ^a) and El^b) that Re and kBT^/Jb increase from their minimum values Re = 
1.2813 and keTj Ji, = 2.679 at A = -2.0 and reach to Re = 1.2932 and keTj Jb = 3.5187 
at A = 0.0 and approaches to its spin-1/2 limit Re = 1.3068 and ksTe/Jb = 5.073 in the 
limit A — )• oo. A similar investigation was made on thin films consisting of a spin-1/2 bulk 
coated with a spin-3/2 surface with using EFT 2J], and it has been reported that Re has 
its maximum value at sufficiently negative crystal fields, and decreases continuously with 
increasing A in positive direction. 

Finally, in Figs. H] and [5], the effect of surface exchange interactions and the presence of 
crystal field interactions on the universality behavior of thin film system have been investi- 
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FIG. 3: Dependence of (a) the special point R^ (b) corresponding bulk transition temperature 
ksTc/Jb as a function of crystal field interaction A. 
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FIG. 4: (Color online) Variation of the shift exponent A for weak surface coupling R = 0.5 cor- 
responding to some selected values of crystal field interaction (a) A = —2.0, (b) A = —1.5, (c) 
A = -1.0 and (d) A = 0.0. 

gated by examining the variation of the shift exponent with A. Only the transition temper- 
atures corresponding to thickness values L {L > 10) which obey the criteria 0.01 < e < 0.1 
are taken into account in the linear fitting process. In Fig. HJ a weak surface coupling such 
as R = 0.5 has been considered. In this case the exponent value is estimated as A = 2.01 
whereas for a moderate ratio of surface to bulk exchange interactions such as R = 1.0, we 
have A = 1.87 which are identical to the values obtained for spin- 1/2 system in Fig. [T](b). 
Hence, we can conclude that the shift exponent A is found to be independent of A. 
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FIG. 5: (Color online) Variation of the shift exponent A for moderate surface coupling R = 1.0 
corresponding to some selected values of crystal field interaction (a) A = —2.0, (b) A = —1.5, (c) 
A = -1.0 and (d) A = 0.0. 

IV. CONCLUSIONS 

In conclusion, with using EFT, we have studied the universal behavior and critical phe- 
nomena in a ferromagnetic thin film described by a spin-1 Blume-Capel Hamiltonian. It 
has been shown that crystal field interaction A plays an important role in layering transi- 
tions. Namely, we have found that the critical value of surface to bulk ratio of exchange 
interactions Re strictly depends on the crystal field interactions A. Strictly speaking, for 
sufficiently negative values of A, Re exhibits a minimum, and continuously increases as A 
increases. In the highly anisotropic limit A — t- oo, it approaches to the value Re = 1.3068 
which corresponds to the critical coupling of spin-1/2 thin film system. In this context, we 
believe that recent calculations reported in the literature should be treated carefully. 

Apart from these, dependence of universality on the surface exchange enhancement, as 
well as the crystal field interactions has been clarified by examining the shift exponent A for a 
wide range of film thickness, 3 < L < 40. In the presence of surface exchange enhancement, 
the exponent A approaches to unity, being independent of crystal fields. In this regard, we 
have concluded that a dimensional crossover may originate as the surface becomes dominant 
against the bulk, and in terms of the exponent A, a ferromagnetic spin-1/2 thin film is in 
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the same universality class with its spin-1 counterpart. 
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